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Abstract 


We analyze the relativistic corrections to annihilation rates of S-wave quarko¬ 
nia within the framework of NRQCD. We show that order v 2 corrections can 
be expressed in terms of the heavy quark pole mass and the quarkonium mass. 
The ratio of hadronic to radiative annihilation rates for % and r] c can there¬ 
fore be predicted accurately. The contributions of color-octet operators to 
the hadronic decay rates of spin-triplet quarkonia are shown to be significant, 
even though they arise at order a 4 in the velocity expansion. We provide a 
rough estimate of the color-octet contributions and extract the value of a s 
from the experimental data on T decays. 
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I. INTRODUCTION 


Since their discovery, heavy quarkonia have been considered an important testing ground 
for quantum chromodynamics [|IJ. By now it is well established that all qualitative features 
of quarkonia (e.g., a confining potential, a positronium-like spectrum, ratios of leptonic to 
hadronic widths) are in agreement with what we expect from QCD. However, in most cases 
we still do not have a fully quantitative description based on first principles. An important 
step towards such a description was made in Ref. . where a formalism of Nonrelativistic 
QCD (NRQCD) was proposed. It is based upon the observation that in a heavy quarkonium 
there are several widely separated momentum scales: the typical kinetic energy of the heavy 
quark, Mv 2 , is much smaller than the inverse size of the quarkonium, Mv, which in turn is 
much smaller than the heavy quark mass M. NRQCD allows one to factor the annihilation 
and production rates for quarkonia into perturbatively calculable short-distance coefficients 
and nonperturbative long-distance matrix elements. This justifies the assumption of “naive 
factorization” for S-wave quarkonium. On the other hand, NRQCD elucidates the role 
of the higher Fock components of the quarkonium wavefunction and explains why naive 
factorization fails for P-wave states 0. 

NRQCD provides a rigorous definition of long-distance matrix elements and thus allows, 
in principle, their calculation on the lattice. Still, given the immaturity of present day 
lattice simulations, one may ask what one can learn from quarkonia without plunging into a 
full-fledged lattice NRQCD computation. What we have in mind here is, first of all, a more 
accurate determination of a s from the ratio of hadronic to electromagnetic widths 0. Such 
a determination in particular could help to clarify the long-standing problem of a possible 
discrepancy between low-energy and high energy measurements of a s [f|]. 

In this paper we analyze the relativistic corrections to the annihilation rates of S-wave 
quarkonia (both spin-singlet and spin-triplet) and apply the results of this analysis to restrict 
the value of a s . The paper is organized as follows. In section || we show that order 
v 2 corrections to the color-singlet part of the annihilation rate can be expressed through 
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the mass of the quarkonium and the heavy quark pole mass. The hadronic widths of the 
spin-triplet states, i/j and T, contain also a piece due to the annihilation of the quark- 
antiquark pair in a color-octet state. In section uni we provide a rough estimate of the latter 
contribution based on the running of the color-octet matrix elements. These results are used 
in section 0 to extract a s from the ratio of hadronic to leptonic widths of T(IS'). 

II. RELATIVISTIC CORRECTIONS TO S-WAVE QUARKONIUM 
ANNIHILATION RATES IN NRQCD 

The Lagrangian of NRQCD p is 

£nRQCD = £ light + £ heavy + SC, (2.1) 

where £n g ht is the usual QCD Lagrangian for gluons and light quarks. £heavy is given by 

£ heavy = (lD t + ^ + X f {lD t ~ X, (2-2) 

with x[) being an operator annihilating a heavy quark, and x being an operator creating a 
heavy antiquark. Both i/s and x belong to the fundamental representation of the color group 
SU(N C ). The last term, SC, includes relativistic corrections to £heavy and is of order v 2 
compared to it. 

The annihilation of the quarkonium is a short distance process (the characteristic mo¬ 
mentum scale is of order M) which, in the framework of NRQCD, is described by adding 
4-quark local operators to SC. The corresponding term has the following structure: 

Vh-,™ = Eilw (2.3) 

n ■*■*■*■ 

where d n is the canonical dimension of O n . The dimensionless coefficients f n ( A) depend on 
the Wilsonian cutoff A needed to define NRQCD and can be calculated by matching the 
NRQCD amplitudes generated by 4-fermion terms with annihilation contributions to the 
scattering in full QCD. Following Ref. [|2|, we take A ~ M. 
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First let us collect the expressions for the decay rates of S-wave quarkonia including the 
first relativistic corrections. According to Ref. the inclusive decay rates of r] c and rj b to 
light hadrons and to two photons are given by 


r(i)c,i - LH) = 2 <^ >4 |o 1 ( 1 s 0 )| I , 0iI ,) 

+ 2 lm J g 4 (1S ° ) MP.eSo)!^} + o(i> 4 r), 
r(i?c,t -> 77) = 2 ^ ( 7 c,t|Oi( 1 S 0 )|>) c , 1 ,) 

+ 2 lm b f 7 4 (lgo> + o(f 4 r), (2.4) 


where Oi( 1 5' 0 ) = 'Pi( 1 Sq) = 1/2 X t (-|D) 2 ^ + h.c 

For the spin-triplet S-states, -0 and T, the situation is more complicated. The leading 
term and the order v 2 relativistic correction to the i + i~ decay rate are proportional to the 
expectation values of C>i( 3 S'i) = ^crx-X^ 0 "^ and V\ ( 3 *S'i) = | —^D) 2 0 + h.c. 

respectively: 


2 Im f ee ( 3 Si) 


r(T tr) = - <t| c?!( 3 5 i_)|t) + — -f;; u (t|^ 1 ( 3 5 , 1 )|t). (2.5) 


2 Im(yf ee ( 3 S'i) 


M 2 


M 4 


The decay rate of 0 or T to light hadrons receives contributions from both color-singlet and 
color-octet components of the quarkonium wavefunction. The color-singlet component can 
only decay into three or more gluons. In contrast, the color-octet component can decay into 
two gluons or into a virtual gluon which then creates a quark-antiquark pair (see Fig. 1.) 
Hence, the color-octet contribution is of order a 2 v A and may compete with the relativistic 
correction to the color-singlet rate which is of order ck 3 u 2 . (We will see in the next section 
that this color-octet contribution is essential for explaining experimental data on T decays.) 
Therefore the inclusive rate to light hadrons is 


T(T -> LH) = 


2 1m / 1 ( 3 5 1 ) 


M 2 

+ r®(T —> LH) 


(T|0 1 ( 3 5 1 )|T) + 


2 Imgifigi) 

M 4 


<T|7M 3 Si)|T> 


( 2 . 6 ) 


The color-octet part of the decay rate r®(Y —> LH) receives contributions from three 
four-quark operators corresponding to the three diagrams in Fig. 1: 
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2 Im (!^P 0 ) + f>M 3 P 2 )) (X | Og( 3 p 0) | X) + 

2 Im^CSo) (X | 0g( i 5o) | X) + 2 lm ^ Sl) (T|0 8 ( 3 5,)|T). (2.7) 

In the latter equation we have used heavy quark spin symmetry to reexpress the expectation 
value of 0 8 ( 3 P 2 ) in terms of O 8 ( 3 P 0 ). The color-octet operators are defined as 

O 8 ( 3 P 0 ) = 1 ^T“ ( -|D • v) xx 'T a (-0 ■ 

OaCS 0 ) =^T a xx'T a 4>, 


r (8) (T -> LH) = 


0 8 ( 3 S 1 )=^aT a x-x^T a ^. 


( 2 . 8 ) 


The short-distance coefficients in the expressions for the decay rates Eqs. (2A-2J\) depend 
on the scheme adopted to define the operators. We wish to derive relations between the 
matrix elements of the operators Pi( 1 5' 0 ), Vi( 3 Si) and Cb^So), (Ti( 3 Si), which can be used 
to express the decay rates including the first relativistic corrections through the leading 
order matrix elements. Since these relations also depend on the choice of the scheme, we 
will discuss this issue in some detail here. We will limit our discussion to the operators 
appearing in the decay rate of the r/b, c - An identical argument can be made for the operators 
appearing in the decay rates of the T and 0. 

The vacuum saturation approximation provides the following estimate for the matrix 
element of Pi( 1 S' 0 ) 


(^7c,&|T > i( 1 5'o)|r7 Ci b) = ^Re (?/ Ci6 |0 t y|O)(O| ((iT>) 2 + X f {iP>) 2 ^\rjc,b) ■ (2.9) 


The operator with two derivatives in this expression is not defined unambiguously beyond 
tree level; for example, one is free to perform a shift 


((0D) 2 X ) f ^ + X f (iD) 2 0) A - (((zD) 2 X ) t 0 + X t(,D) 2 0) A (2.10) 

+C(A, M) (xV) A - 


Here A is the Wilsonian cutoff, and C(A,M ) is a power series in a s starting with a term 
of order a s . (In what follows we will omit the subscript A, with the understanding that all 
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operators are regularized using the Wilsonian cutoff.) Among all possible definitions of the 
operator with two derivatives, only a subset satisfies the NRQCD velocity counting rules. 
According to these rules, the matrix element of the operator with two derivatives should 
scale as p 2 relative to the matrix element of yW as p —> 0: 


<0| ((^D) 2 y) t V' + y t (fD ) 2 ^(p, —p) 

(o|x¥Ip, -p) 


o(p 2 ). 


( 2 . 11 ) 


Imposing this condition removes the freedom to redefine the operator as in Eq. ( ]2 .1 C| ) . 

With the convention Eq. ( p.ll|) it is particularly simple to determine the short-distance 
coefficient of the operator Oi( 1 S'o) at next-to-leading order (NLO) by comparing the v -1 and 
v° contributions in full QCD and NRQCD. On the NRQCD side it is sufficient to calculate 
the annihilation of a quark-antiquark pair via the operator 0 1 ( 1 S'o). The choice Eq. fl2.11|) 
guarantees that there are no NLO contributions from the operator Vi^Sq) proportional to 


v 1 or v°. 


Other ways of defining the operators are, of course, possible. If the p —> 0 limit of 
the expression on the right-hand side of Eq. (|2.11| ) is nonzero, the operator Ri( 1 S'o) will 
contribute to the annihilation rate at order v and v°. In this case one needs to know 
the leading order short-distance coefficient of Vi^Sq) in order to determine that of 0\ (’Ao) 
at NLO. In fact, unless one requires operators with arbitrarily many derivatives to satisfy 
conditions similar to Eq. (|2.11 ), they all have to be taken into account in the computation 
of the short-distance coefficient of C>i( 1 S'o) to next-to-leading order. 

In Ref. J2[] it was assumed that the operators with two or more derivatives need not be 
taken into account when computing the NLO short-distance coefficient of OiQSo). In other 
words, it is implicit in Ref. || that the scaling behavior of the operators is given by the 
NRQCD counting rules or, equivalently, that Eq. ( 2.1 1|) and similar equations for operators 
with more derivatives are satisfied. We adopt the same convention here. 

The equations of motion for the quark fields to leading order in v 2 are 


iDt + 


D 2 ' 
2 M 


■0 = 0 , 


iD, 


D 2N 
2 M 


X = o. 


( 2 . 12 ) 
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They can be used to trade the spatial derivatives in the operator with two derivatives in 
Eq. (|2.9|) for time derivatives acting on the quark fields: 


((®) 2 x) V + X* (®) 2 V> + j4(A, M) x ] ii> = 2 M id, (xty) • 


(2.13) 


Here A is a scheme dependent coefficient whose expansion in powers of a s starts, in general, 
with the term of order a s . The term proportional to x'i> has to be included in the relation 
Eq. ( |2.13|) , because x'V ; mixes into the operator with two derivatives under shifts as in 
Eq. (|2.10|) . There are corrections to Eq. ( |2.13|) at higher orders in v 2 , but for our purposes 
it is sufficient to take into account only the terms shown. 

Let us show that with the convention Eq. ( [2,1 lp A(A,M) is zero to all orders in a s . 
Evaluating Eq. (|2.13p between vacuum and a quark-antiquark state |p, —p) yields: 


((?.D) 2 y) t ^ + y t (*D) 2 ^|p,-p) = (2p 2 - H(A, M)) (0|xty|p,-p)- (2.14) 


We have used the identity 

(o| id t (xH) Ip, -p) = <0| [xV> h 


Ip.-p) = ^(o|xVIp,-p), 


(2.15) 


where H is the NRQCD Hamiltonian. The identity Eq. (|2.15|) holds to all orders in a s 
because the asymptotic state |p, —p) is an eigenstate of H with eigenvalue p 2 /M. Dividing 
both sides of Eq. ( |2.14| ) by (0|xV|P, — p), taking the limit p —> 0 and using Eq. ( |2.11|) , one 
sees that H(A, M) = 0. 

Using Eq. ( |2.13|) with A = 0 yields the following relation: 

{Vc,b\'Pi( 1 So)\r]c,b) = ^Re(?7c ) 6|V’ t x|0)(0| ((iD) 2 xf ip + X t (iD) 2 ijj\r] c>b ) 

= MRe <?7 c ,6| / 0 t x|O) (01 id t (xV) I Vc,b) (l + 0(u 2 )) 

= ME^irj^O.CSo)^) (l + 0(v 2 )) . (2.16) 


Here E Vcb is the energy of the quarkonium state. Since in NRQCD the rest mass of the 
quarks is not included in the energy of a quarkonium state, we can express E Vc b in terms of 
the mass of the rjb tC and the quark pole mass 
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(2.17) 


Er]c,b — ^»?6,C 2M, 


Eqs. (12.16| ) and (p . 1 7|) and similar relations for the operators in the spin-triplet decay rates 
allow us to express the relativistic corrections to the decay rates in terms of the leading order 
matrix elements and the ’’binding energy” of the quarkonium M qq — 2 M q . For example, the 
decay rate of the T to light hadrons now takes the form 


T(T -> LH) = 


2 Im/ipS*, 


M 2 

+ T (8) (T -> LH), 


(T|0 1 ( 3 S' 1 )|T) 1 + 


M r -2Mlmg 1 ( 3 S 1 )' 
M Im / 1 ( 3 S' 1 ) 


(2.18) 


and similar expressions hold for the other decay rates Eqs. ( pi.4|) and Eq. 

All coefficients in Eqs. ( p.4|) and Eqs. (|2.5H2.7[ ), except / 8 ( 3 P 0 ), Zs( 3 P 2 ), and gi( 3 Si), have 
been calculated to the necessary order in a s in Ref. 0. The coefficients / 8 ( 3 Po), ZsP-^) can 
be extracted from Eqs. (A9-A13) of Ref. 0: 


3n(N 2 -4) 


n(N 2 - 4) 


Im/ 8 ( 3 Po) = ^f^^(M), Im/ 8 ( 3 P 2 ) = v f AT - a 2 s (M). 


(2.19) 


AN C sy n 5N C 

To extract 5 'i( 3 S'i) to leading order in a s we need to compute the three-gluon annihilation 
rate of a free quark-antiquark pair to order v 2 in their relative velocity. Fortunately, this 
computation has already been performed in the context of e + e _ annihilation ||. We have 
checked the results quoted in these papers. The annihilation rate of a free quark-antiquark 
pair in a spin triplet state to order v 2 turns out to be 

0 197r 2 - 132 


T (ggpS'i) -> 3 g,vj = T (gg( 3 ^i) -> 3^,0 


+ 0(v 4 


( 2 . 20 ) 


12vr 2 - 108 

where T (qq( 3 Si) —> 3 g,v) represents the annihilation rate of the quark-antiquark pair in a 
spin-triplet state, with v denoting the velocity of the quark in the center of mass frame. 
Comparing with the corresponding amplitude in NRQCD, we obtain 

ImtqpS'i) 197T 2 — 132 


(1 + 0(a s )) 


( 2 . 21 ) 


Im/iPSi) 12 tt 2 -108 
One consequence of the last equation is that for the spin-triplet states the order v 2 
relativistic correction to the hadronic rate is unexpectedly large. For the 6-quark pole mass 

















in the range 4.6 — 4.9 GeV, the correction to the T(IS') decay rate can be as large as 25%, 
and still bigger for radially excited states. For m^ ole ~ 1.3 GeV, the correction to the J/-0 
hadronic decay rate is about 150%. Its magnitude makes one question the usefulness of the 
nonrelativistic expansion for charmonium. For spin-singlet states, the relativistic corrections 
are of the expected size. 


III. ESTIMATES OF THE COLOR-OCTET MATRIX ELEMENTS 


In order to use our expressions for phenomenological applications estimates for the color- 
octet contributions are needed. Following Ref. 0, we can obtain very rough estimates by 
solving the renormalization group equations for the color-octet operators. To order v 4 and 
leading order in a s we find 


a4_(T|O s (‘S„)|T) 

a4-( T |0 8 ( 3 S,)|T) 

aA(T|0 8 ( 3 P„)|T> 


0, 

4(Af, 


A)a s 


N c ttM 2 
4 Cfcx 


81 N c tt 


(T|0 8 ( 3 Po)|T>, 


(M r — 2M) 2 (Y|C> 1 ( 3 S' 1 )|T), 


(3.1) 


where 


ifij^lT) = M 2 (M t - 2M) 2 (T|0,( 3 S 1 )|T> 


(3,2) 


was used in Eqs. O- We can express the matrix elements at the factorization scale A ~ M 
in terms of those at a low scale A ~ Aq CD by solving Eqs. (|3.1|) . The color-octet operators 
mix between themselves and with color-singlet operators. Formally, the terms coming from 
the mixing with color-singlets are logarithmically enhanced. To get a rough estimate of the 

color-octet matrix elements, we assume that these terms dominate. This yields: 

/ \ \ 2 


(T| 08 ( 3 Sl )| T) » 8(iV V~ 4 i Cf (Mt ~ . 2M)2 


(T|Os( 3 P„)|T) 


81V^7T 2 M 2 

(Mr — 2M) 


3d)IT\~ 4 Cf (1 \/r_. 0/1^2 2?r ln 


A) \ol s (M ) 

1 


81A5.7T 


A) \ a s(M), 


JJ (T|C> 1 ( 3 S' 1 )|T), (3.3) 

(T|C? 1 ( 3 <S' 1 )|T). (3.4) 
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In the same spirit we set (Y|C> 8 ( 1 So)|T) ~ 0 since it does not acquire a logarithmically 
enhanced contribution. 

In order to check whether these estimates are reasonable, we consider the following “ratio 
of ratios”: 


T (T(mS') -► LH) /r (T(mS') Plr) 
~ T (T (nS) -»■ LH) /r (T (nS) -»• FF) ' 


Substituting Eq. ( p. 7|) and a similar expression for the dileptonic rate into Eq. (|3.5|) we get 


M m — M n r® (Y(mS) —> LH) 
M m r (T(mS) -> LH) 


(T(raS) -> LH) 
T (T(nS') -*• LH) 


+ 0(r> 4 , ayr? 2 ). (3.6) 


Neglecting the color-octet contribution completely, we obtain f ? 12 — 1.5, in disagreement 
with the experimental value = 0.95 ± 0.15 [| 6 ! ]. In order to evaluate the color-octet 
contribution numerically, we need the pole mass of the b quark. Various methods give M in 
the range 4.6 — 4.9 GeV [|7[|8| corresponding to the theoretical value Ru — 0.95 — 1.53. For 
the lower quark masses our estimates shift the value of f ? 12 much closer to the experimental 
number. We take this to be an indication that our estimates give reasonable order of 
magnitude values for the color-octet matrix elements. 


IV. APPLICATION TO THE DETERMINATION OF 

As we have seen in the previous section, the leading relativistic correction to the anni¬ 
hilation rates of r] c and r /5 is expressed in terms of the 5-quark pole mass. The latter can be 
extracted from the measurement of moments of the photon spectrum in the inclusive b —> sy 
decay [9j, from inclusive semileptonic b —► c decays [7j, or from sum rules for quarkonia [Q . 
Therefore the ratios of hadronic to radiative decay rates of r/ c and rib are ideal for determin¬ 
ing a s . Unfortunately, these measurements are very hard (though not impossible) to do. [] 

Where is still one problem on the theoretical side. Next-to-leading order (NLO) perturbative 
corrections to these ratios are very large |2],[l(J, and one would like to know NNLO corrections to 
have some idea about the convergence of the perturbation series. 
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The Ts are much easier to study from the experimental point of view, but the theoretical 
interpretation is complicated by the presence of the color-octet contribution. We use the 
results of the previous section to estimate this contribution to the hadronic decay rate. For 
M = 4.6 — 4.9 GeV the color-octet contribution ranges from 0 to 9%. Therefore we take 
9% as an estimate of the error from neglecting it. The order v 2 relativistic correction is 
evaluated with M in the same range as above. We use the renormalization scale dependence 
of the NLO prediction to estimate the error from perturbative NNLO corrections. 

Having adopted such estimates of theoretical uncertainties, we use the experimental value 
T (Y(IS') -> LH) /T (T(1S) -»• g+l~) = 37.3 ± 1.0 @, to determine a s (M) = 0.154 - 0.218. 
This corresponds to 


a s (M z ) = 0.097- 0.117 


(4.1) 


at the scale M z . The higher value of a(M z ) corresponds to the lower value of the 5-quark 
pole mass. This range for a s (M z ) overlaps with the 1 a confidence interval of the LEP 
measurement a s (M z ) = 0.120 ± 0.004 [TT[] . The accuracy of our extraction being limited 
by theoretical uncertainties, the range in Eq. ( |4 . 1| ) should not be interpreted as a 1 a error. 
We do not quote here the values of a s obtained from T(2S') decays because the theoretical 
uncertainties are much larger, and also because the accuracy of data on T(2S) is worse. 

Further improvements in this determination of a s would come from a more accurate 
extraction of the 5-quark pole mass, and also from a NNLO perturbative calculation of the 
short-distance coefficient Im/ 1 ( 3 S'i). For example, knowledge of the pole mass to within 
50 MeV would reduce the uncertainties roughly by a factor of two. 

We have shown that order v 2 relativistic corrections to annihilation rates of the S-wave 
quarkonia can be expressed in terms of the quarkonium “binding energy”. For spin-singlet 
states this observation makes it possible to predict accurately the ratio of hadronic to ra¬ 
diative decay rates in terms of a s and the heavy quark pole mass. However, a calculation of 
the NNLO perturbative contributions to the short-distance coefficients is necessary to en¬ 
sure that perturbative corrections are under control. For spin-triplet states, which are much 
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more accessible experimentally, the color-octet component of the quarkonium wavefunction 
may contribute signihcantly to the hadronic annihilation rate, although the corresponding 
contributions are of order u 4 in the nonrelativistic expansion. Therefore, knowledge of the 
expectation values of color-octet operators is needed, if we want to predict the hadronic to 
leptonic ratio for spin-triplet quarkonia. We used crude estimates based on renormalization 
group equations to deduce the uncertainties due to the color-octet contributions. From the 
experimental data on the Y(IS') decays we extract a s (Mz) = 0.097 — 0.117, the major part 
of the uncertainty coming from the uncertainty in the 6-quark pole mass. Further experi¬ 
mental and theoretical efforts are needed to obtain a better estimate of a s from quarkonia 
decays. 
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FIG. 1. Contributions of higher Fock states to the hadronic annihilation of spin-1 S-wave 
quarkonia. All diagrams shown here contribute to the rate at order 
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